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ABSTRACT 

We investigate the effect of ram-pressure from the intracluster medium on the strip¬ 
ping of HI gas in galaxies in a massive, relaxed, X-ray bright, galaxy cluster at z = 0.2 
from the Blind Ultra Deep HI Environmental Survey (BUDHIES). We use cosmological 
simulations, and velocity vs. position phase-space diagrams to infer the orbital histo¬ 
ries of the cluster galaxies. In particular, we embed a simple analytical description of 
ram-pressure stripping in the simulations to identify the regions in phase-space where 
galaxies are more likely to have been sufficiently stripped of their HI gas to fall below 
the detection limit of our survey. We find a striking agreement between the model pre¬ 
dictions and the observed location of Hi-detected and non-detected blue (late-type) 
galaxies in phase-space, strongly implying that ram-pressure plays a key role in the 
gas removal from galaxies, and that this can happen during their first infall into the 
cluster. However, we also find a significant number of gas-poor, red (early-type) galax¬ 
ies in the infall region of the cluster that cannot easily be explained with our model 
of ram-pressure stripping alone. We discuss different possible additional mechanisms 
that could be at play, including the pre-processing of galaxies in their previous envi¬ 
ronment. Our results are strengthened by the distribution of galaxy colours (optical 
and UV) in phase-space, that suggests that after a (gas-rich) field galaxy falls into 
the cluster, it will lose its gas via ram-pressure stripping, and as it settles into the 
cluster, its star formation will decay until it is completely quenched. Finally, this work 
demonstrates the utility of phase-space diagrams to analyze the physical processes 
driving the evolution of cluster galaxies, in particular HI gas stripping. 

Key words: Galaxies:clusters:general - Galaxies:clusters:invidivual (Abell 963) - 
Galaxies: evolution 


1 INTRODUCTION 

Much work has been done to understand the physical 
processes driving galaxy evolution. It has been shown by 
many observational and theoretical studies that the envi- 


galaxies, (e.g. Dressier 19801: Lewis et al. 20021: Gomez et al. 

20031: KooDmann & Kennev 20041: Boselli & Gavazzil 20061: 

Poggianti et al. 

2006:; 

Moran et al. 20071; Bamford et al. 

2009; Jaffe et al. 

2011b 

). An unassailable result from ob- 


servations of galaxies in the local universe is that cluster 
galaxies are more likely to have an elliptical or SO mor¬ 
phology than field galaxies, which are typically star-forming 


spirals dDresslejIl98d ). In addition, it is known that the 
morphology of galaxies correlates with their stellar popula¬ 
tions. In a colour-magnitude diagram, elliptical and SO (i.e. 
early-type) gala xies follow a tight correlation (the so-calle d 
“red sequence”; iBau ml ll959l : IVisvanathan fc Sa ndagell 1 9771). 
that i ndicates they have old stellar populations I Bower et al.l 


ll992l ; lArag6n-Salamanca et al JfT 99 3l : I .Jaffe et al, 1201 lal ). Spi- 

ral and irregular (i.e. late-type) galaxies instead tend to dis¬ 
play a wide range of bluer colours, often referred to as the 
“blue cloud”. The existence of a clear colour bimodality in 
clusters strongly supports the idea that environment plays 
an important role in quenching the star formation in galax¬ 
ies. Moreover, observations of distant galaxies have shown 
a steady rise in the fraction of SOs in clusters, along with 
a comparable decline in the spiral fraction since z ~ 0.5, 
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while the fraction of elliptical galaxi es remains constant 
llFasano et al.l l2000l : iDesai et al .II 2 OO 7 I) . When considering 
the predicted evolution o f large-scale str ucture from cos¬ 
mological simulations (e.g. iDe Lucia et al.ll20ll) . this result 
suggests that spiral galaxies transform into SOs as they tran¬ 
sition from the field into clusters. The physical mechanism(s) 
driving such a transformation is however still under debate. 


A particularly important ingredient for understanding 
galaxy evolution is the atomic gas. Neutral hydrogen (HI) 
not only provides a reservoir from which H 2 forms, fueling 
star formation, but is also a sensitive tracer of different en¬ 
vironmental proc esses, such as r am-pr essure stripping and 
tidal interactions. iGunn fc Gotti dl972l) first predicted (ana¬ 
lytically) that ram-pressure can be effective in removing the 
galaxies’ interstellar medium (ISM) as they pass through 
a dense intracluster medium (ICM), thereby quenchin ; 


Bravo-Alfaro et al. 

to 

o 

o 

£2- 

20011; Poggianti & van Gorkom 

2001 

: Kennev et al. 

2004|; 

Crowl et al. 20051: 

Chung et al. 

2007 

, 2009|; Abramson et al. 201 ll: Scott et al 

2010. 2012!: 


is indeed disturbed and eventually truncated and e xhausted 
in clusters. Recent simulations (see [Roedigerti2009l . for a re¬ 
view) have further supported the idea that this happens 
via r am-pressure stripping and gravitational interactions 
(e.g. IVollmer et al.ll200ll:lyollmeiil2003l : lTonnesen fe Brvanl 
120091 : Kapferer et al. ~ l2*009h . It is important to note that 
ram-pressure is expected to be more effective at the cen¬ 
tre of massive clusters, because of the higher densities 
of the ICM, and the high velocities of the galaxies. In 
fact ram-pressure stripping models indicate that HI can 
be fully depleted after a transit through the ICM’s core 
dRoediger fc Bruggenll2007l 'l. Moreover, the galaxies in mas¬ 
sive clusters also suffer fast encounters with other clus¬ 
ter members. T he accumulat i on of such encounters, known 
as harassment dMoore et al.l H 99 J ) . can cause disk thick¬ 
ening and gas fuelling of t he central region . Harassment 
can strip HI from galaxies dDuc fc Bournaudl 120081 ) and is 
most effective in low-surface-bright ness gala xie s, an d to¬ 
wards the core of galaxy clusters dMoore et al.l 119991 9. Fi¬ 
nally, it has also been suggested that galaxies can be pre- 
proc essed in less massive group s before they fall into clusters 
(e.g. IZabludoff fc Mulchaev|l99ji. T his idea has recen tly re¬ 
ceived observational dB alogh et al.lll999l: iLewis et al . 2002 : 


IVerdes-Montene gro et al.||200l|; Jaffe_et_ajj j20ljj_ and theo- 

retical support (jMcGee et alJl200fll : lDe Lucia et al . 2M). In 

order to acquire a complete understanding of the physical 
processes affecting galaxies through their lifetimes, it is thus 
necessary to couple their observed properties to the cosmic 
environments which they have inhabited throughout their 
lifetimes. 


Due to technological limitations, studies of the HI con¬ 
tent of galaxies have so far mostly been carried out in the 
local universe (< 0.08). To investigate the relation between 
gas reservoirs and star formation, as well as its evolution, 
we carried out a Blind Ultra Deep HI Environmenta l Sur¬ 
vey IBUDHIES lVerheiien et aLll2007i : I.Taffe et al.ll2013l ). that 
focuses on two pencil-beam pointings containing 2 galaxy 
clusters at 2 ~ 0.2, and the large-scale structure around 


them. The clusters, Abell 2192 and Abell 963 (A2192_l and 
A963-1 from now o r[]i are very different from each other. 

In Ijaffe et ahl d 20121) we studied A2192_l, and found 
that this intermediate-ma ss cluster has a sig nificant amount 
of substructure (see also Ijaffe et~aP 120131 ) . The incidence 
of Hl-detections significantly correlates with environment: 
at large clustercentric radii (>2-3 virial radii), where the in¬ 
falling substructures were found, many galaxies are detected 
in HI, while at the core of the forming cluster, none of the 
galaxies are Hi-detected. We found that this effect starts to 
become significant in low-mass groups that pre-process the 
galaxies before they enter the cluster, so that by the time 
the group galaxies fall into the cluster they may already be 
HI deficient or even devoid of HI. 

A963_l, on the other hand, is a massive cluster with 
significant X-ray emission centred on the brightest cluster 
galaxy (BCG). This cluster is thus ideal to study the inter¬ 
action between galaxies and the ICM. In this paper, we focus 
on A963_l to study the effect of ram-pressure in the strip¬ 
ping of HI gas from cluster galaxies, and the later quenching 
of their star formation. 

Rec ent works by Mahajan^t all ( 201 lh : Oman_et i jil 


(2013); IHernandez-Fernandez et al.l (l2014h : iMuzzin et al 


( 2014 ); iTaranu et al.l 1 20141) have highlighted the useful¬ 
ness of phase-space diagnostics to infer the assembly his¬ 
tories of clusters. In simulations it is possible to trace back 
the orbital histories of galaxies in clusters, and associate 
different populations (e.g. infalling, backsplash, or galax¬ 
ies in the virialized part of the cluster) to different phase- 
space locations. Simulations can provide a 6-dimensional 
phase-space (all 3 components of both the position and 
velocity vectors: x,y,z,v x ,v v ,v z respectively). Observations 
of distant galaxies, however, only give 3 projected param¬ 
eters: 2 coordinates, and the line-of-sight velocity (a, 5, 
and vios respectively). For simplicity we refer to the sim¬ 
ulated phase-space as “3D”, and to the observed one as 
“projected”. Despite the loss of information, it is still pos¬ 
sible to infer the assembly history of cl uster galaxies from 
their projected phase-space (PPS ; see e.g. lOman et al.ll2013l : 
IHernandez-Fernandez et all 20141 ). The uncertainties intro¬ 
duced by the projection effects can be minimized by utilizing 
large (mass-limited) samples of galaxies in regular clusters of 
similar mass and redshift. Moreover, this approach certainly 
provides a more complete view of environment compared to 
the standard ID view (galaxy properties as a function of 
radial clustercentric distance only). 

In this paper, we investigate the environmentally driven 
evolution of galaxies in A963_l, a massive cluster at z = 0.2, 
making use of the unique BUDHIES dataset. We exploit the 
Hl-data and ancillary data to map the gas stripping and 
later suppression of the star formation activity of galaxies 
in phase-space. To interpret our data, we compare it to a 
simple ram-pressure model built into a cosmological simula¬ 
tions, that predicts the regions in phase-space where HI gas 
is expected to be stripped from the galaxies in our cluster. 

The paper is organized as follows: In Section Owe sum- 


1 As in Jaf fe et all (120131) , we adopt the names “A963_l” and 
“A2192_l” to distinguish these clusters from the other structures 
we found along the same line-of-sight, such as galaxy groups and 
cosmic sheets. 
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marize the data collected for BUDHIES fSec. l2.il) . the main 
characteristics of the cluster fSec. 12.21) . and the properties of 
its member galaxies fSec. l2.3ll . In Section[3]we show how the 
cluster galaxies are distributed in PPS, highlighting the lo¬ 
cation of the Hl-detections. In Section[XT]we create a simple 
ram-pressure stripping model, which we embed in the cos¬ 
mological simulations presented in Section [472] These are, in 
turn, used to derive the assembly histories of the observed 
clusters and predict the final location of stripped galaxies in 
PPS. In Section l4~3l we combine our observed and simulated 
results and discuss the effect of ram-pressure stripping, and 
compare our results with previous work. We also use galaxy 
colours in the PPS analysis as tracers of post-stripping evo¬ 
lution in Section m to better understand the late stages of 
star-formation quenching as a function of the orbital histo¬ 
ries of the cluster galaxies. We summarize our findings and 
draw conclusions in Section [5] 

Throughout this paper we assume a concordance 
ACDM cosmology with =0.3, Qa =0.7, and Ho =70 
km s _1 Mpc -1 . 

2 OBSERVATIONS 

2.1 BUDHIES data 

BUDHIES is an ultra-deep HI survey of galaxies in and 
around two clusters at 0.16 ^ 2 ^ 0.22. The survey has 
an effective volume depth of 328 Mpc and a coverage on 
the sky of ~12xl2 Mpc for each cluster, which allows us to 
study the large-scale structure around the clusters. At the 
core of our study are the ultra-deep HI observations, carried 
out at the Westerbo r k Synthesis Radio teles cope (WSRT 
lYerheiien et al.l 120071 : lYerheiien fe et all l2010h from which 
we have detected 127 galaxies in A963_l and 36 in A2192_l 
with HI masses > 2 x 10 9 Mq. In addition to the HI data, we 
have obtained B- and R-band imaging with the Isaac New¬ 
ton Telescope (INT), NUV and FUV imaging with GALEX 
(Fernandez et al. in preparation), Spitzer imaging, and near- 
infrared imaging with UKIRT. SDSS photometry is avail¬ 
able for the galaxies in the held, and in a limited number 
of cases, SDSS spectra are also available. Finally, we have 
carried out a spectroscopic campaign with WIYN and the 
AF2 on the William Herschel Telescope (WHT), from which 
we have assessed cluster memberships and characterized the 
environment in the surveyed volumes (1 Jaffe et al.ll201 j i. 

2.2 A963_l 

With a dynamical mass of 1.1 x 10 15 Ji -1 Mq (<t c / = 993 km 
s -1 ), A963_l is the most massive cluster in BUDHIES. It 
is found at z = 0.2039 and its R 200 was estimated to be 
1.55 h -1 Mpc djaffe et al .11 201 jh . The cluster’s total X-ray 
luminosity is Lx — 3.4 ± 1 x 10 44 ergs -1 /i -2 dAllen et al.1 
l2003i) . Its regular and centrally concentrated X- ray emissio n 
is cen tered on a central cD galaxy (see Fig. 6 of ISmith et alJ 
21 )orj;, . suggesting A963_l is a relaxed cluster. However, in 
J affe et al.l (1201311 we perfo rmed a Dressler-Shectman test 
1 Dressier fc Shectman|[l98a 'l on the spectroscopic members 
of A963_l out to larger radii, and concluded that the “re¬ 
laxed” part of the cluster is at its core, while there is a mild 
presence of substructure in the outskirts. In addition, small 


X-ray subs tructures have bee n observed at r > R 200 (see 
Fig. B.6. in lZhang et al.ll2007l . see also Haines et al. in prep, 
for a full description of the X-ray substructures). In this pa¬ 
per, we focus on the global properties of A963_l, ignoring 
small inhomogeneities in the ICM. This is justified because 
the vast majority of the X-ray flux is centrally concentrated 
and very regular. The substructure in the cluster outskirts 
does not dominate. We leave the detailed study of the HI 
properties of galaxies as a function of X-ray structure for a 
future paper. 


2.3 The cluster galaxies 

In the following, we present derived properties of the mem¬ 
ber galaxies of A963_l. 


2.3.1 Stellar masses 

To compute stellar masses (M*) for all the spectroscopically 
confirmed cluster members, we first converted the 5-band 
SDSS photometry available to absolute magnitudes using 
our spec troscopic red shifts, and the interpolation utility In¬ 
ter Rest ll Taylor et al .1 120091 ). We then computed th e stella r 
masses following the prescription bv IZibetti et al.l ([2009) , 
with a Kroupa initial mass function. 

The distribution of M* for the cluster members is shown 
in Figure[l] In the figure, log(M*) is plotted against R mag¬ 
nitude. Hi-detected galaxies are marked with a blue square 
and occupy the region of lower stellar masses and fainter 
R magnitudes. The horizontal dashed red line indicates the 
stellar mass cut (M* = 10 10 Mq) used to define a mass- 
limited sample, adopted in some of the analysis presented 
in this paper. This mass cut roughly matches the magni¬ 
tude cut where our spectroscopic ca mpaign is most c om- 
plete (dashed orange vertical line, see I Jaffe et ali (l2013l l for 
details on the target selection). 


2.3.2 Colours 

To distinguish different galaxy populations, we constructed 
colour-magnitude diagrams (CMDs) using the UV and op¬ 
tical data. I 11 the left panel of Figure [2] we show the B — R 
vs. R CMD for A963_l. We fitted the red sequence using the 
reddest (B — R > 1.9) early-type galaxies (i.e. elliptical and 
SOs). The best ht (B — R = —0.042 x R + 2.98) is shown 
by the solid black line in the plot. We further define the 
blue cloud to contain all galaxies bluer than 0.25mag below 
the red sequence (blue filled symbols below dashed line). We 
note that the morphological classification used to define the 
red sequence was the product of visual inspection of the INT 
B and R images. Several co-authors assigned morphologies 
to the galaxies and these were later combined. The quality 
of the images, however, was too poor to robustly assign mor¬ 
phologies to all the cluster galaxies. In fact we could only 
determine a rough morphological type for a small sample 
of galaxies. To avoid introducing scatter to our analysis, or 
having to reduce the sample significantly, in this paper we 
only used the morphological information to ht the red se¬ 
quence. In the analysis however, we use colours instead to 
broadly separate star-forming, late-type galaxies (blue) from 
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A963_1 




Figure 2. Left: B-R vs R diagram of galaxies in in A963_l. The red sequence fit is shown by the solid line. Red-sequence galaxies (red 
diamonds) are defined to be at most di0.25mag bluer than the red sequence fit (dotted line). The blue-cloud galaxies (blue circles), are 
the ones that lie below the dotted line. Right: NUV-optical colour-magnitude diagram showing the red-sequence and blue-cloud galaxies 
defined in the left panel, and dotted lines separating regions (colour ranges) that correspond to different stellar-population ages, as 
labeled (see text for description). In both panels HI detections are highlighted with open blue squares. 



Figure 1. R-band magnitude versus log(M*) for the cluster 
members (black circles). The Hl-detections are distinguished with 
open blue squares. At the top and right-hand side of the plot, 
histograms of the R-band magnitude and log(M*) are also shown. 
The dashed red line indicates the mass cut adopted in the analysis 
of Section 14.31 (see also Section and the dotted blue line is 
the mass of the Milky Way-like galaxy used as a prototype to 
explore the effects of ram-pressure stripping in the cluster (c.f. 
Section l4.ll) . The magnitude limit adopted in the spectroscopic 
campaign is indicated by the vertical orange line. 


those with quenched star formation and early-type morphol- 
ogy (red). 

The colour classification shown in the left panel of Fi¬ 
gure [2] provides an indication of the ages of the stellar popu¬ 
lations in these galaxies. The majority of red-sequence galax¬ 
ies, for instance, are expected to have stopped forming stars, 
while blue galaxies are still hosting star formation and thus 
contain a significant fraction of young stellar populations. 
In the Figure [5] we have also highlighted the location of 
Hi-detected galaxies. As expected, these populate the blue 
cloud. 

Tabic [l] lists the sizes of the total “blue” and “red” 



No. gals. 

No. massive gals. 

No. Hl-detec. 

Red 

128 

109 

1 

Blue 

118 

49 

65 


Table 1. The galaxy subsamples used in this paper, split by their 
optical colours, as defined in the left panel of Figurc[2] The colums 
are: sample name (“Red” and “Blue”, corresponding to the red 
sequence and blue cloud), number of galaxies, number of galaxies 
in the mass-limited sample (M* > 10 10 Mq), total number of 
Hi-detected galaxies. 

galaxy samples, as well as the number of galaxies in the 
mass-limited and Hi-detected samples. 

The right panel of Figure[3]shows the UV-optical colour- 
magnitude diagram. UV is very sensitive to the ages of young 
stellar populations. It is therefore helpful to utilize NUV- 
optical colours as a tool to roughly infer the ages of the 
young stellar component in the galaxies. We do this by in¬ 
specting the NUV-r colour evolution in a single-stellar po¬ 
pulati on model of solar metallicity from lBruzual fc Charlotl 
(|2003i rl. We adopt the same x-axis range and symbol colour- 
code as in the left panel to easily locate the red-sequence and 
blue-cloud galaxies previously defined in optical bands. As 
expected, the galaxies spread over a larger range in colour in 
the UV-optical colour-magnitude relations. We choose three 
NUV-r colour bins that correspond to the ages of a single 
stellar population of < 300 Myr; 300—600 Myr; > 600 Myr, 
at the redshift of our target cluster. These bins are identi¬ 
fied by the dotted lines in the right panel of Figure 0 We 
emphasize that these lines should only be treated as good 
reference values, and not exact ages of the stellar popula¬ 
tions. We note that we did not consider internal extinction 
in our computation, but we use these colour cuts only as a 
lower limit for a given age of the young single stellar po¬ 
pulation. We will use these colour bins in the analysis of 
Section roi 

2 We note that o ur results do not cha nge if we use the models 
from Starburst99 (iLeitherer et al.ll 19991) 
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3 A PHASE-SPACE VIEW OF A963_l 


To investigate the dynamical state of the cluster under study 
and infer (statistically) the orbital histories of the differ¬ 
ent galaxy populations, we have constructed line-of-sight 
velocity vs. radial position phase-space diagrams for the 
spectroscopically-confirmed members, as shown in Figured] 
The x-axis of the plot shows the projected distance from the 
cluster centre, r, normalized by R200, while the y-axis dis¬ 
plays the peculiar line-of-sight velocity of each galaxy with 
respect to the cluster recessional velocity (An), normalized 
by the velocity dispersion of the cluster. In particular, we 
plot the absolute value, taking advantage of the symmetry 
in the velocity distribution of the cluster. This is: 

|Au| = \c{z-z d ) | 

a (1 + Zd)cr 

where z is the redshift of a given galaxy, z c i the redshift of 
the cluster, and c the speed of light. 

One can think of the position of galaxies in PPS as lower 
limits to their real 3D radius and velocity values. In other 
words, a galaxy cannot be at a smaller radius as it appears 
in projection. Likewise, it cannot have a smaller velocity 
difference than its line-of-sight velocity difference. 

The grey contours in Figure [3] show the galaxy number 
density distribution in phase-space. This plot represents a 
snapshot (at the observed redshift) of the cluster assembly 
process, where infalling galaxies come towards the cluster at 
large radii, gaining velocity as they reach the cluster centre. 
They will then settle into the cluster, partly by means of dy¬ 
namical friction (which is especially important for high mass 
galaxies), creating the dense “wedge” in PPS, as we will see 
in Section H21 The overall distribution of cluster galaxies in 
PPS thus arises from “virialization”, and it is constrained 
by the radial dependence of the escape velocity in PPS. This 
is shown by the dashed black line, th at corre sponds to the 
escape velocity in a lNavarro, Frenk fc Whitel (119961 . NFW) 
halo, projected to the line-of-sight. This is: 


Where 


V5SP r < r ™ 

otherwise 


2GM200 

3i?200 s 


( 2 ) 


I< = g c 


ln(l + cs) 


- ln(l + c) + 1 


(3) 


s = 


-R200 


7r r 
2 R200 


(4) 


we have further normalized v to the velocity dispersion of 
A963_l. Galaxies inside the dashed black line in Figure [3] are 
expected to be gravitationally bound to the cluster. 

The “virialized” region of the clu ster is roughly locate d 
at r < R200 and Av < 1 . 5 a (see e.g. IMahaian et al.ll201l| ). 
although the shape of the virialized region is best approxi¬ 
mated by a triangle, rather than a square (as will be shown in 
Figures |4] [5] [7] an dEJ. 70% of the galaxies inside this region 
have been inside the cluster (r^n < R200) for more than 4 
Gyrs, which corresponds to the crossing timcQof our cluster. 
In other words, the virialized region contains galaxies that 
have most likely experienced many pericentric passages. 

In the PPS of Figure El we have distinguished low-mass 
galaxies from more massive ones with symbol sizes, and 
we can see that low- and high-mass galaxies are roughly 
evenly distributed across PPS in the cluster (although there 
are regions notably dominated by massive galaxies, such as 
0-3 < r/J ?200 < 1 and |Au|/<r < 0.7). This is not neces¬ 
sarily a real effect but perhaps reflects the spectroscopic in¬ 
completeness of the low-mass galaxy sample. Although we 
always plot all galaxies for reference, it is important to be 
aware that a significant assessment of the incidence of differ¬ 
ent galaxy populations can only be done in the mass-limited 
(or spectroscopically complete) sample. For this reason we 
mostly focus on the high-mass galaxies, plotted with bigger 
symbols in our analysis. 

In Figure El we also highlight, for the moment, HI- 
detected galaxies (blue squares), as these observations are 
the core of the survey. Later, in Figures [7] and [9] we will 
see how other galaxy properties are distributed in PPS. It is 
clear that the Hi-detected galaxies are located outside the 
cluster virialized region and that there is also a notable lack 
of Hl-detections outside the virialized region at high veloci¬ 
ties and small radii. We indicate with the solid purple line in 
Figure El the boundary between the “Hi-rich” region (right 
of the line) and the region devoid of Hl-detections (or “HI- 
poor”; left of the line). Note that this line was drawn by 
eye, as reference. The cause of this segregation between the 
Hi-detected and non-detected galaxies will be thoroughly 
investigated in the following sections. 

We emphasize that all Hi-detected galaxies were plotted 
in this diagram, no HI mass cut was assumed. However, 
the effect of the primary beam attenuation of the WSRT 
translates into a variation of the Hi-mass limit with distance 
from the pointing centre. We will describe these effects in 
Verheijen et al. (in preparation). It is remarkable that we 
see strong decline of the incidence of Hi-detected galaxies 
towards the cluster centre (see also Figure El despite this 
effect, which makes HI harder to detect at larger radii. 


Sc = 


ln(l + c) 


-1 


c 


(1 +c) 


(5) 


and c is the concentration parameter, which we fixed 
at c = 6 (the exact value of c does not change the results 
significantly). Note that the 3D radius scales with the pro¬ 
jected r, with scatter, as: r = cos (tp)r 3D , where ip is the 
projection angle, that can vary between 0 and n/2. There¬ 
fore, on average, r ~ §^30 ■ We also projected the 3D escape 
velocity ( v esc ) to a line of sight escape velocity (u) using the 
average conversion factor. Inevitably, the projections (i.e. in 
both v and r) lead to increased scatter in PPS. I 11 the figure, 


4 THE IMPACT OF RAM PRESSURE 
STRIPPING IN PPS 

4.1 Ram-pressure stripping model 

To explain the segregation of Hi-detected galaxies in PPS 
(Section El, we explore the effect of ram-pressure stripping 
in A963_l. We do this by creating a simple model, based on 


3 A “crossing time” refers to the time that an infalling galaxy 
takes to cross from r = R200 to the other side of the cluster. 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 

r /^200 


Figure 3. The observed PPS for all galaxies in A963-1. Galaxies in the mass-limited sample (M* > 1 () 10 A/q) are represented with 
darker and larger filled symbols, Hi-detected galaxies are enclosed by a blue open square, and the grey countours follow the number 
density of the galaxies. The dashed black line corresponds to the escape velocity in a NFW halo (see text for details). The solid purple 
line delimits the region (right of the line) containing the vast majority of the Hl-detections in the cluster. 


the analytic description prov ided in ,£rjmn ! ^_£ro£tj ((lgZjl). In 
particular, following IHernandez-Fernandez et al.l (|2014l ). we 
compute the “intensity” (rf) of ram-pressure (Pram) exerted 
by the intra-cluster medium on an infalling galaxy, as the ra¬ 
tio of such pressure over the anchoring self-gravity provided 
by the galaxy (n goi ). That is, r/ = Pram/P-gai- 

We model the ICM gas density profile of the cluster wit h 
the standard /3-model ( Cavaliere fe Fusco-Femianol 1 19761 ). 


Picm ( r 3 d ) — Po 


1 + 


R c 


1 - 3 / 3/2 


(6) 


A963 l’s ICM model 

P 

0.5 

R c 

65 kpc 

PO 

14.7X10 -3 cm -3 


Table 2. Parameters of the /3-model for A963_l, used to compute 
the intensity of ram-pressure as a functio n of distance from the 
cluster centre. The values were taken from lRizza et al . 1 Gill). 

where r' is the radial distance from the centre of the 
galaxy, R d is the disk scalength, and So the central surface 
density of the galaxy: 


where, as before, the 3D radius scales with the projected r, 
with scatter, as: r ~ -r, 

7 7T 3D 

Now, the ram-pressure is defined as: 

Pram = PlCM V gal ( 7 ) 

Again, v ga i is the velocity of the galaxy in 3D, which 
we also project to the line-of-sight velocity, v. 

We then model the galaxy’s anchoring pressure, Il ga i, 
as the following: 


II g a i — 27rGS a E g (8) 

Where E s and S g are the density profiles of the stellar 
and gaseous disks respectively. We assume an exponential 
profile for the disk: 


S = X 0 e- r ' /Rd 


(9) 


So = 


M d 

2 *Rl 


( 10 ) 


where M d is the mass in the disk. _ 

We adopt the /3-model parameters from Iftizza et al.l 
(119981) to compute P ra m for A963_l. These values are sum¬ 
marized in Table [2] 

Finally, we construct a disk galaxy model to test the 
effect on ram-pressure on it. The values used to model the 
galaxy’s anchoring force (listed in Table 0 are a reason¬ 
able representation of a Milky Way (MW)-like disk galaxy. 
Throughout the rest of this paper we thus refer to the model 
as a “MW” model. Note that the stellar mass of this test 
galaxy corresponds to the most typical mass in our mass- 
limited sample (see dotted blue lines in Figure [lj. Using 
these values we obtain Yl ga i = 1.24 x 10 _11 N m -2 at the 
centre of the galaxy (i.e. at r' = 0, where the restoring force 
is the stron gest). N ote that this value i s sma ller than the 
one used bv IHernandez-Fernandez et al.l (120141 ). 
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value 

referece 

A^d, stellar 

4.6x IO 1U M 0 

Bow & Rix (2013) 

A/d,gas 

~0.06 X Md t stellar 

Gavazzi et al. (2008) 

^d, stars 

2.15 kpc 

Bow & Rix (2013) 

Rd , gas 

1.7 X.Rd,stars 

Cavatte et al. (1994) 


Table 3. Values used to compute the restoring force in the disk 
of a MW-type galaxy. 


We can now predict the region in PPS where ram 
presure is expected to strip the gas off the MW-like galaxy in 
A963_l, i.e. the case where rj = 1. Jumping ahead to the bot¬ 
tom panel of Figure[4](or Figures[5][T] and[9j, we can see the 
case of complete stripping (r' = 0) plotted as a grey dashed 
line in the PPSs of Figures 14151 [3 and [9] Galaxies infalling 
into the cluster (from the right-hand side of the line towards 
the cluster centre), will eventually cross this line, losing their 
gas reservoirs. This line thus constitutes a good measure of 
the presence of HI in PPS. However, because our HI obser¬ 
vations have a detection limit of ~ 2 x 10 9 Mg, there will 
be some galaxies that will fail to be detected before the gas 
disk is truncated all the way down to r' = 0. In other words, 
we expect to stop seeing Hl-detections at larger r/I ?200 than 
the grey dashed line predicts. To account for this effect, we 
have re-computed Ii ga i at the truncation radius at which a 
Milky Way galaxy would fall below BUDHIES’ HI detection 
limit. The truncation radius, rt > 0, is computed from the 
remaining gas mass in the galaxy, / = ———, following: 


/ 


1 + 


e ~ rt/Rd 



( 11 ) 


Solving for rt we obtain a value of 6.1 kpc, which yields 
the new H ga i at r' = rt, shown by the dashed green line in 
Figure |T] (and Figures [5] [7] and 0. Galaxies to the left of 
this line are thus expected to have been stripped enough to 
fall below our detection limit. We call this region the ram- 
pressure stripping cone. 

For future reference and application to other systems, 
in Appendix m we have explored how much the cluster and 
galaxy model affect the mapping of the effect of ram-pressure 
stripping in PPS. Interestingly, we find that the effectiveness 
of ram-pressure in PPS does not vary much with galaxy mass 
(except for the most massive galaxies). We will come back 
to this in Section m 

We also note there are several caveats to our ram- 
pressure stripping model, that need to be considered for 
the analysis. To start, our model depends on at least 6 
variables (3 spatial and 3 in velocity), and we are probing 
its effects in two dimensions. Some scatter in PPS due to 
projection effects (in both radius and velocity) is thus ex¬ 
pected. We also note that, despite the simplicity of Gunn 
& Gott’s prescription, their model has been validated by 
many simulations. Although it is well-known that their for¬ 
mula is best at predicting the stripping radius in the case of 
face-on stripping, other inclination angles (with respect to 
the i nfall trajectory) c a n lead to more com plex behaviour 
fe.g. lAbadi et al.l Il999l : IVollmer et al.ll200ll ~). However, the 
amount of stripping is o nly signi ficantly al tered for nearly 
edge-on galaxies (iRoediger fe Bruggenll2006h . Finally, obser¬ 
vations and simulations of galaxies in clusters have revealed 
that it is possible to have gas stripping at larger radii than 


Example galaxy orbits in phase-space 



Figure 4. The orbit tracks in 3D (top) and projected (bottom) 
phase-space of two massive galaxies randomly extracted from the 
simulations. The points correspond to different time-steps in the 
simulation. The first point in each case is labeled with tj, which 
corresponds to 7.942 Gyr before the final snapshot, tf. In the 
bottom panel, the dashed grey line indicates the area in PPS 
(left of the line) where a MW galaxy is expected to be completely 
stripped as it falls into the cluster. The dashed green line instead 
delimits the area where the model galaxies are stripped enough 
to fall out of the detection limit of our survey. The dashed black 
line indicates the region in PPS where galaxies are most likely 
to have been in the cluster for more than a pericentric passage 
(i.e. the “virialized” region). The region to the right of all dashed 
lines thus contains galaxies that have most likely recently joined 
the cluster. Clearly, at ti Galaxy A was falling in to the cluster 
for the first time, while Galaxy B had aready crossed i? 200 - As 
time passes, galaxy B oscillates in position and velocity drawing 
a “wedge” in phase-space (inside the “virialized” region). This is 
what will happen to galaxy A after enough time has passed. 


what Gunn fe G otti 1 1972|j predict fe.g. iKennev et aDl2004 
iTonnesen et al.l 12007 ). due to the presence of non-smooth 
substructure in the ICM, which translates into scatter in 
ram-pressure strength at fixed radius. As mentioned in Sec¬ 
tion E3 we will combine the BUDHIES data with deep X- 
ray observations to thoroughly investigate this scenario in a 
future paper. 

Despite all the aforementioned caveats, in Section 14.31 
we will show how the distribution of Hi-detected galaxies in 
PPS is in exceptionally good agreement with the prediction 
from our model. 


4.2 Cosmological simulations 

The recipe described in Section 14.11 is a good measure of 
the velocity and radial dependence of stripping for galaxies 
on their first infall. However, in PPS we have a mixture 
of galaxies that are infalling, and others that have been in 
the cluster for several pericentric passages. To predict the 
location in PPS where we expect to find stripped galaxies 
(virialized or infalling), we apply our simple ram-pressure 
model to dark matter halos orbiting in the potential of a 
cluster formed in a cosmological model. The advantage of 
doing this is that one can trace back the galaxies to study 
how they build up the cluster (and thus PPS). 

We use the high resolution cosmological N-body 




























Projected 



Figure 5. The 3D (left) and projected (right) phase-space of the simulated cluster galaxies. The top panels show galaxies that are 
infalling into the cluster for the first time. They clearly occupy the region in phase-space close to the cluster “envelope” and far away 
from its core. The middle panels show all (virialized and infalling) galaxies, which occupy the entire trumped-shaped area of phase-space. 
At the bottom, we show only massive galaxies (i.e. with masses similar to that of our MW prototype). In the right-hand panels we 
further indicate some regions of interest: the ram-pressure stripping cone where galaxies will fall below our HI detection limit (left of the 
green dashed line), the cone of total HI stripping (left of grey dashed line), and the region in PPS where most galaxies are expected to 
have been in the cluster for more than a pericentric passage, i.e. the “virialized” region. 


simulations, fu l ly de scribed in IWarnick Sz Knebel ((2006) ; 
[Warnick et al.l (|2008h . which we re run with the a daptive 
mesh refinement code MLAPM (jKnebe et al.ll200lh . These 
simulations used approximately one million dark matter par¬ 
ticles, each with a mass of 1.6 x 10 8 fl _1 Mg, while the spa¬ 
tial force resolution was ~ 1 kpc. The halos were found 

us ing AHF, a par allelized version of the MHF algorithm 
of lGill et al.l (120041 '). The minimum particle number allowed 
for a halo was set at 20, leading to a minimum resolved halo 
mass of - 3 x 10 9 Mq. 

To emulate A963_l as much as possible, we cho se the 
cluster “Cl” from Table 1 of lWarnic.k fe Knebel (j2006f ). as it 
is the most massive cluster in these simulations. The virial 
mass and radius for that cluster are: AI v i r =2.9 x 10 14 M o , 
and R v i r = 1.355 Mpc, and the formation redshift zj = 
1.052 Gyr, which is equivalent to an age of ~ 8 Gyr. 

We follow the orbits of individual halos in the cluster, 
as shown in the examples in Figure 0] This figure not only 
shows the typical trajectories of galaxies in phase-space as 
they fall (and eventually settle) into the cluster, but also 
reveals the distortions that can appear due to projection 
effects. For each halo, we measure the time-evolving ram- 
pressures, assuming the /3-model presented in Section 14.11 


(see Table [2] in particular). We assume the halos contain a 
MW-like disk with the same properties listed in Tablc[3] and 
flag them as “stripped” (i.e. non-detected in HI) once they 
have been truncated down to rt and have thus fallen below 
the detection limit of our survey. 

In the top panels of Figure [5] we plot the final positions 
(at z=0) in a phase-space diagram, using only the halos that 
are infalling into the cluster for the first time, and distin¬ 
guish the ‘HI detected’ (open light green circles) from the 
‘non-detected’ or stripped galaxies (filled dark green circles). 
It is noticeable how galaxies drop out of our detection limit 
as they pass through the ram-pressure stripping cone de¬ 
scribed in Section rm (area left of the green dashed line 
in the right-hand panels of Figure EJ. This effect is visi¬ 
ble in both 3D and projected phase-space, although it is 
less distinct in projection. In the middle panels of the figure 
we plot all halos (infalling and virialized), and see that the 
non-detected galaxies are not only in this cone, but elon¬ 
gate along a “wedge”, as indicated by the dashed-black line 
in PPS. This line is not an exact boundary, but a rough indi¬ 
cation of the region where galaxies have a higher probability 
of having been in the cluster for more than a pericentric pas¬ 
sage. 
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r /-^200 


Figure 6. Distribution in PPS of the probabilities that a MW- 
type galaxy has been stripped enough to be undetected by 
BUDHIES in a cluster with the ICM model of A963_l. The red- 
orange areas trace the regions where galaxies have most likely 
been stripped of their gas. The thick black line is the same as 
in Figure [3] that shows the boundary between the region where 
most Hi-detected galaxies are (right of the line) and the region 
without Hl-detections (left of the line). 

From Figure Owe can see how the position of the simu¬ 
lated galaxies in PPS evolve with cosmic time. On first infall, 
the galaxies distribute along the escape-velocity caustic we 
saw in Figure O avoiding low r and low v values (top panels 
of Figure 0). As the cluster evolves and becomes virialized, 
the galaxies cluster in the low r and low v region, as shown 
in the middle panels of the Figure (filled darker symbols, see 
also example orbits in Figure [dj- This shows the capability 
of PPS to predict, in a statistical sense, the time since infall 
of cluster galaxies. 

In order to ensure that dynamical friction is of the 
correct magnitude in our results, we subsample our ha¬ 
los, creating a ‘MW-like’ subsample with initial halo mass: 
0.25 x 10 11 — 7.5 x 10 12 Mq. This approximately selects ha¬ 
los with the correct mass to match the stell ar mass of our 
MW-like model (see Eq. 3 in lGuo et aI.ll20fol ). There are 96 
in total in this mass range (bottom panel of Figure [5]). By 
subsampling our halos, we are also assuring that there is a 
good match between the disk mass and the halos, although 
the general trends observed in the middle and bottom panels 
of Figure [5] are not dissimilar. To be safe, however, we only 
use the sample plotted in the bottom panels in our analysis. 

To project the 3D phase-space, we randomize the line- 
of-sight to the cluster, spinning the cluster by a random 
angle about the x-axis, then another random angle about 
the y-axis. The line-of-sight is then chosen down the z-axis, 
with r = \/x 2 + y 2 , and vi os = v z . 

Then, we compute the probability, per PPS bin, that 
a galaxy has been stripped. To do this, we spin the cluster 
1000 times and build up an average in each pixel of the PPS 
diagram (fraction of HI non-detected to the total number 
that lie in that pixel). This is shown in Figure E] where red- 
orange regions trace the highest probabilities for a galaxy 
to have been stripped. In the figure, we have overplotted 


the line separating Hl-detections from non-detections in the 
observed PPS shown in Figure [3] The model predictions for 
where the effect of ram-pressure is strongest in PPS is in 
striking agreement with our observations, qualitatively re¬ 
producing the region in PPS where there are no Hi-detected 
galaxies (left of the line). Although the results shown in Fi¬ 
gure [5] were computed specifically for MW-type galaxies in a 
massive cluster (like A963_l), and assumed the HI detection 
threshold of BUDHIES, this PPS probability map can be 
used in other cluster studies as a guideline for where ram- 
pressure will be at play, provided the galaxy and cluster 
properties are not dramatically dissimilar from this study 
(see Appendix O- The virialized “wedge” will not change 
significantly. What will mostly change is the location of the 
ram-pressure stripping cone, that depends on the truncation 
radius (or HI detection limit). 

4.3 The effect of ram pressure in HI stripping 

In this section, we combine the results of the previous sec¬ 
tions, to link the effect of ram-pressure stripping with the 
orbital histories of galaxies in the studied cluster. 

Figure [7] shows the PPS of galaxies in A963_l. As in 
Figure [3] galaxies of different masses are shown with dif¬ 
ferent symbol sizes, where the larger symbols correspond to 
the mass-limited sample. Hl-detections are highlighted as 
before. We also distinguish between red-sequence and blue- 
cloud galaxies (defined in Figure [2|. As in previous figures, 
the area to the left of the dashed green line indicates the 
ram-pressure stripping cone where galaxies will fall below 
our detection limit (see Section f4.ll) . while the grey dashed 
line delimits the region of complete gas stripping. The area 
below the black dashed line instead, represents the “virial¬ 
ized” zone, and the area outside the stripping and “virial¬ 
ized” regions thus contains galaxies that have recently joined 
the cluster, as labeled. In Table[4] we further list the fraction 
of different galaxy populations (for the mass limited sample 
only) in all the regions of PPS defined. 

The most striking result of Figure [7] is that both the 
stripping and virialized areas are almost completely devoid 
of Hl-detections, and that their combined areas remarkably 
coincide with the region where simulated galaxies that have 
been stripped of their gas are located in PPS (left of line 
in Figure [6]), regardless of their time since infall. Note that 
the stripped and “virialized” regions are not only devoid of 
Hl-detections, but are also dominated by red galaxies (espe¬ 
cially when considering the “virialized” region without the 
overlapping stripped cone, see Table II. This is consistent 
with galaxies aging as they spend more time in the cluster. 
The “Recent infall” region on the other hand is full of blue 
galaxies, likely coming from the field, as well as a significant 
fraction of red cluster members. 

Taking the mass-limited sample, we find that only a mi¬ 
nority of the few blue galaxies in the “stripped” and “virial¬ 
ized” regions (or “Hi-poor” region) are detected in HI (9%). 
Inversely, outside these regions (in the “recent infalls” zone), 
we see that the majority of the blue cluster galaxies are de¬ 
tected in HI (67%). Our results strongly suggest that field 
galaxies (typically blue, gas rich, star-forming spirals) that 
enter a massive cluster, will experience HI gas stripping on 
their first infall (i.e. first passage through the ram-pressure 
stripping cone), and that the intensity of the ram-pressure 
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Figure 7. The same PPS as in Figure 171 is shown, but with added information: As before, galaxies in the mass-limited sample (M* > 
1O 1O M0) are represented with larger filled circles, Hi-detected galaxies are enclosed by a blue open square, and grey countours follow 
the number density of the galaxies. In addition, galaxies are divided into red-sequence (red diamonds) and blue cloud (blue circles), as 
defined in the top panels of Figure [2] The dashed grey line delimits the area in PPS (left of the line) where a MW galaxy is expected to 
be completely stripped as it falls into the cluster. The dashed green line delimits the area (to the left of the line) where the model galaxies 
are stripped enough to fall out of the detection limit of our survey (see Section |4.2D . As in previous figures, the additional dashed black 
line indicates the region in PPS where galaxies are most likely to have been in the cluster for more than a pericentric passage (i.e. the 
“virialized” region). The region to the right of all dashed lines thus contains galaxies that have most likely recently joined the cluster. 



Figure 8. Fraction of Hi-detected galaxies in A963_l as a func¬ 
tion of r/i ?200 (solid blue line). The blue square at (arbitrarily) 
large radii corresponds to the field fraction in BUDHIES, for com¬ 
parison. 


they experience is sufficient to make them undetectable by 
BUDHIES. 

Although our galaxy sample is very incomplete at M* < 
10 10 Mq it is interesting to note that the “HI-poor” region 
seems to be the same for low and high mass galaxies. This 
is consistent with the results shown in Appendix El where 


we found that variations in the mass of the model galaxy 
does not change significantly the predicted regions in PPS 
where ram-pressure stripping is most effective (except for 
very massive galaxies that are harder to strip). These two 
findings thus suggest that the orbital histories of the galax¬ 
ies in the cluster are more important than their mass, to 
determine the amount of stripping they have experienced. 

As mentioned above, there is a noticeable fraction of red 
(and HI-poor) galaxies in the “recent infall” region of PPS 
(59%, see Table 2J- as well as some blue galaxies that are 
not detected in HI (33%). The presence of stripped galax¬ 
ies in the cluster outskirts is not a new finding. In the Virgo 
cluster, st rongly stripped HI disks have been observed in the 
outskirts (Crowl & Kenney |2004 . Their morphologies and 
stellar populations, together with an observational estimate 
of their orbital histories, suggest these galaxies have been 
strip ped in the outskirts of clusters, before crossing the cen¬ 

tre JChung et all2009l : ICrowl fe Kennevl2008l : iKennev et al.l 
12004 . 

Using cosmological simulatio ns of cluster formation and 
evolution, iTonnesen et al.l d2007ll found that although gas 
stripping occurs primarily in the central region of clusters, 
it is an important mechanism out to the virial radius. This is 
due to the wide scatter in ram-pressure strength that galax¬ 
ies can experience at fixed radius. They also point out that 
the timescale for complete gas removal is ^ 1 Gyr, and that 
galaxies in the cluster periphery (r > 2.4 Mpc) often accrete 
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cool gas. More recent simulations bv lCen et al.l J2014h fur¬ 
ther highlight the fact that the cluster environment can have 
an impact on the gas content of galaxies out to large radii, 
and that galaxies lose their cold gas as they fall into clusters, 
within a single radial round-trip. It is still unclear, however, 
whether the observational evidence for ram-pressure strip¬ 
ping outside clusters is related to the non-smooth density 
distribution of the ICM. 

An alternative explanation for the presence of quenched 
galaxiefl in the cluster outskirts (in PPS) is that not all 
of the galaxies in the outskirts of the cluster are necessar¬ 
ily falling in for the first time, but backsplashing instead. 
The amplitude of the backsplash after first infall and the 
time they will take “virializing” depend on a number of 
factors, which makes the identification of backsplash galax¬ 
ies in PPS very difficult. In fact, the regions where these 
galaxies tend to be found also host other galaxy popula¬ 
tions, such as galaxies o n t heir f irst infall or older cluster 
members (e.g. lOman et al.112013 1. This means that back- 
splashing is not a dominant effect at any particular PPS 
location. In particular, backsplashing is expect ed to be_in- 
signif i cant at large radii (r > 2 x R 200 , see e.g. iBahe et al.l 
120131:1 Muriel fe Coendall201.fi . so it cannot account for all 
of the observed quenched fraction in the “Recent infall” re¬ 
gion. Moreover, our own simulations, did not show signi¬ 
ficant quenching in this area of PPS (Figure [6J , when we 
applied the ram-pressure stripping prescription to it. This 
implies that not all red galaxies outside the “virialized” re¬ 
gion can be attributed to backsplashing. An additional effect 
is needed to account for the high fraction of red galaxies out 
to the largest r. One likely possibility is that some galaxies 
had already been quenched by the time they fell into the 
cluster (e.g. though group pre-processing). T his idea is sup¬ 
ported by the findings of lHaines et al.l (|2007f ). who showed 
(using SDSS) that the fraction of passively-evolving galax¬ 
ies in the field can be as high as 50%, when considering 
the most luminous systems. It is thus possible that some of 
these red galaxies in the outskirts of A963_l were already 
red when accreted into the cluster from the field or a lower- 
density environment. I 11 fact, it is possible that a portion 
of these galaxies could have been “pre-processed” (e.g. via 
mergers) in smaller groups that were accreted into the clus¬ 
ter. This idea is supported by our previous study of the 
distribution o f Hi- d etecte d galaxies in the forming cluster 
A2192_l (i.Iaffe et al.ll2012l b We will investigate this further 
in future papers by studying the Hi-content of galaxies in 
the small infalling groups visible in X-rays (see Section ^. 21) . 

We also com pare our results with those found by 
IBahe et al.l (I2013T) . who recently analyzed the (hot and cold) 
gas content and star-formation rate of galaxies as a func¬ 
tion of environment with hydrodynamic cosmological simu¬ 
lations. They find a significant decrease in the HI content 
and star formation with (3D) distance from cluster centre 
(r 3D ) out to very large distances (r 3D ~ 5 x R 200 , which 
is ~ 3 x R 200 in projection, with scatter). Leaving aside 
the phase-space view for a moment, and only concentrating 
on r/R 2 oo, we get the trends sho wn in Figure El that can 
be compared with those found by IBahe et al.1 (l2013l b The 


4 By “quenched galaxies”, we are referring to those galaxies with 
quenched star-formation 


blue line corresponds to the fraction of Hi-detected galaxies 
corrected for HI incompleteness caused by the WSRT’s pri¬ 
mary beam attenuation. The correction was co mputed using 
the HI mass function of iMartin et al.l (l2010l l as reference. 
We note how ever t hat there might be a variation between 
IMartin et al.l (l2010h ’s field HI mass function at low redshift, 
and that of A963_l. We ignore this effect when calculating 
the corrections for Figure[8] as we expect the difference to be 
within the plotted Poisson errors. An in-depth assessment 
of the HI mass function in BUDHIES is beyond the scope 
of this paper. We plan, however, to investigate the HI mass 
function as a function of both redshift and environment in 
a future paper. 

Despite the mentioned caveat, Figure [8] shows that the 
fraction of Hi-detected galaxies in the cluster clearly in¬ 
creases with r/R. 2 oo, until it reaches a plateau that matches 
our measurement for the field, within the errors. The gen¬ 
eral trend seen in BUD HIES is in agreement with those pre¬ 
dicted by IBahe et al.l J2013h (the red line in their Figure 2 
corresponds to a cluster mass similar to that of A963_l), 
when considering that their r 3D range (0 < r 3D < 5 x R 200 ) 
shrinks to our plotted range in projection. Although our HI 
fraction reaches the field value at smaller r, we note that the 
comparison can only be done qualitatively, as the absolute 
values of the fractions depend on their specific recipes and 
assumed thresholds. Measurements of Hl-deficiencies in clus¬ 
ter galaxies at low z also sho w a significan t environmental 
effect only out to ~ 2 R 200 llSolanes et al.ll200lh . in agree¬ 
ment with our results. Ou r resu lts ho wever are not neces¬ 
sarily in disagreement with IBahe et al.l (l2013ll ’s conclusions, 
that evoke a combination of ram-pressure stripping, back¬ 
splashing and group-preprocessing to explain the observed 
trends. 

4.4 Star-formation quenching post HI stripping 

Even if field galaxies are affected by ram-pressure stripping 
as they first fall into the cluster, as proposed in 14.11 an in¬ 
falling galaxy that crosses the BUDHIES ram-pressure cone 
(green dashed line) could continue to form stars in the re¬ 
maining gas disk. Only when the galaxy has been completely 
stripped by the ICM (if it crosses the dashed grey line, pass¬ 
ing very near the core of the ICM, or if it crosses less dense 
regions of the cluster sufficient times), will its star formation 
rate rapidly decline, moving the galaxy from the blue cloud 
to the red sequence. At the same time, the galaxy will move 
in PPS until it settles into the virialized area. This scenario 
is supported by the properties of the blue galaxies inside the 
stripping cone, that have redder (although still blue) optical 
colours and on average, an earlier-type morphology (from 
the few available morphologies in the sample) with respect 
to the blue galaxies in the “recent infall” region. Moreover, 
when stacking the spectra of the few blue galaxies in the 
stripped zone for which we have AF2/WYFFOS data, we 
see clear signs of H$ absorption, which suggests declining 
star-formation activity. 

Turning our attention to the star-formation in the 
stripped and non-stripped galaxies, we now analyze the loca¬ 
tion of UV-detected galaxies in PPS. In particular, it is very 
interesting to see where the bluer galaxies (in NUV — r ) are 
located in PPS with respect to the rest. This is shown in Fi¬ 
gure [9] (large circles). Similarly to the Hi-detected galaxies, 
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Figure 9. The same PPS of Figure[7]is shown, but only for massive galaxies (M* > 1O 1O M0), this time colour-coded (and differentiated 
by symbols) by the stellar-population ages, inferred from their NUV — r colours (see Section 12.3.2b , as labeled. Hl-detections are 
highlighted with blue open squares. Lines are as in Figure JT] 
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Table 4. Fraction (in %) of different galaxy populations in the different regions of PPS. Following the column order: (1) the ram-pressure 
stripping cone (left of green dashed line in figures), (2) the virialized region (below black dashed line), (3) the combination of both of these 
regions (which coincides with the “Hi-poor” region of Figure [3}> (4) the ram-pressure stripping cone region excluding the overlapping 
area of the virialized zone, and (5) the region mostly populated with galaxies recently accreted to the clusters (right of dashed lines, 
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refers to the population of Hl-detections within the blue galaxies. The last 3 rows show: fuv hiue ? the fraction of galaxies in the bluest 
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reddest bin, that also contains galaxies that are not detected in NUV (> 600 Myr), as defined in Figure [2] All fractions were corrected 
for spectroscopic completeness, considering the mass-limited sample only. 


the galaxies with the bluest NUV — r colours (filled purple 
circles) tend to be located outside of the stripping cone, and 
are scarce in the “virialized region” (see Table Q}. This is 
consistent with the findings of Section 14.31 since the blue¬ 
ness of these galaxies suggests that they host a very young 
stellar component (< 300Myr), as discussed in Section [2.3.21 
When we consider slightly redder galaxies (those with a stel¬ 
lar population component that is ~ 300 —600Myrs old, filled 
turquoise squares), we see that they occupy a larger area of 
PPS: although they are common outside the stripped area, 
they can also be found in the ram-pressure stripping cone, 
at the higher velocities, suggesting they have not been in the 
cluster for a long time. It is remarkable that the very-young 
and young galaxies are less common in the “virialized” area 


of PPS (see Figure [9] and Table [4]). Finally, if we consider 
the redder (older) galaxies (smaller symbols), it is clear that 
the distribution in PPS is less segregated, although they are 
predominant in the virialized region. The vast majority of 
these galaxies (including the ones not detected in NUV) are 
red-sequence galaxies (see Figure [5J, so it is not surprising 
that their distribution in PPS mimics that of red galaxies 
in Figure [7] As mentioned earlier, it is possible that some 
of these galaxies are indeed falling in to the cluster for the 
first time, but were quenched (or started to be quenched) in 
their previous environment. 

In summary, the distribution of UV-detections in PPS 
proved to be a good tool to study star-formation as a func¬ 
tion of the galaxies’ orbital histories. The results presented 
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in this section, together with those presented in Section \4 . 3 1 
support a scenario in which blue, star-forming, late-type 
galaxies falling into the cluster from the field, are stripped 
of their HI gas (down to at least BUDHIES’ detection limit) 
on their first passage through the ICM. After the gas is en¬ 
tirely removed, the star-formation is reduced and eventually 
quenched. This causes the “virialized” part of the cluster 
(in PPS) to be populated by red and “dead” systems. In ad¬ 
dition, our analysis clearly shows that not all galaxies take 
this evolutionary route. In fact, the presence of a significant 
fraction of gas-poor galaxies with old stellar populations in 
the infall region of the cluster requires an additional mech¬ 
anism other than ram-pressure stripping (see discussion in 
Section roi) . 


5 SUMMARY AND CONCLUSIONS 

We use BUDHIES to investigate the role and timescale of 
HI gas stripping and later decline of the star-formation of 
galaxies in A963_l, a massive regular cluster at 2 = 0.2. We 
do this by constructing line-of-sight velocity vs. projected 
position phase-space diagrams of the cluster galaxies out to 
the cluster outskirts (~3x R 200 )- Using cosmological simu¬ 
lations, these diagrams are then used to infer (statistically) 
the orbital histories of the cluster galaxies. 

Having a basic understanding of the cluster’s assem¬ 
bly history allows to more effectively study the evolution of 
the cluster galaxies themselves, in particular, the transfor¬ 
mation (s) they may experience as a consequence of having 
joined the harsh cluster environment. We investigate the ef¬ 
fect of environment on the galaxies by studying the phase- 
space distribution of their HI, optical and NUV properties. 
The presence/absence of HI is an excellent indicator of envi¬ 
ronmental effects that lead to gas removal (e.g. ram pressure 
stripping), whilst the optical and UV colours can be used to 
distinguish galaxies with recent star-formation from those 
with old stellar populations. 

To interpret the observed segregation of Hi-detected 
and non-detected galaxies in PPS, as well as the distribu¬ 
tion of blue/red galaxies, we model th e effect of ram- p ressur e 
stripping using the the prescription of I Gunn fc Gotti (1 19T2I) . 
and adopting the HI mass limit of our survey (~ 2 x 10 9 Mq), 
the ICM density profile of A963_l, and the stellar masses of 
the observed cluster galaxies. We embed this model into the 
cosmological simulations so that we can study the effect of 
ram-pressure stripping as a function of the orbital histories 
of cluster galaxies. We then predict the regions in phase- 
space where galaxies are expected to be stripped of their 
gas. In particular, we create a map that shows, per phase- 
space bin, the probability that a galaxy has been sufficiently 
stripped to be undetected by our survey. The predicted map 
coincides strikingly well with the location of BUDHIES’ HI- 
detected and non-detected blue (late-type) galaxies. 

Our findings strongly suggest that HI gas in late-type 
galaxies can be stripped (down to at least the detection 
limit of our survey) during their first infall into the cluster. 
In phase-space, the stripping will occur once they have ap¬ 
proached the dense ICM core and/or gained enough velocity 
to cross the “stripping” area in phase-space (small cluster¬ 
centric distances and/or high line-of-sight velocities). After 
this, they will oscillate in phase-space (as they bounce back 


and forth in the cluster, gaining and losing speed) until they 
settle near the cluster centre due to dynamical friction. This 
long (> 4Gyrs) process gives rise to the dense “virialized” 
region observed in PPS. During this time, the galaxies’ gas 
reservoirs are expected to be completely exhausted, which 
in turn causes their star formation rates to be quenched. 

To test this scenario, and to study the evolution of the 
stellar populations of galaxies after gas stripping, we com¬ 
bine the Hl-observations with optical (B — R) and NUV 
(NUV — r) colours in our analysis. Although we find that 
galaxy colours also segregate somewhat in phase-space, the 
trends are less clear than when using HI. However, the “viri¬ 
alized” part of the cluster is dominated by red-sequence 
galaxies, as expected, while the cluster outskirts or infall re¬ 
gion has a mixed population of galaxies. Moreover, we find 
that optically blue galaxies are less rare than Hi-detected 
galaxies in the “stripped” region of phase-space. This is con¬ 
sistent with blue, star-forming galaxies first losing their gas 
reservoirs via ram-pressure stripping, and slowly declining 
their star-formation activity, as described above. In this sce¬ 
nario, the “stripped” galaxies (i.e. the ones not detected by 
BUDHIES) will still have blue colours for a while, until the 
gas is completely removed. Once this happens, the gas-poor 
galaxies will transit to the red sequence. The distribution 
of the galaxies’ NUV colours in phase-space further support 
this scenario. In particular, we see how the bluest galax¬ 
ies in NUV — r (i.e. the ones with younger stellar popu¬ 
lations) mimic the distribution of the Hi-detected galaxies, 
while galaxies with “green” UV colours (intermediate-age 
stellar populations) are not only present in the cluster out¬ 
skirts but are also found in the “stripping” region of phase- 
space. Finally, the galaxies with the oldest stellar popula¬ 
tions (reddest NUV colours or undetected in NUV) domi¬ 
nate the “virialized” region of the cluster. 

The described pathway to quench the star-formation in 
galaxies is, however, not unique. In fact, an additional in¬ 
sightful result from our analysis is the finding of a significant 
population of passive (red, gas-poor, early-type) galaxies in 
the outskirts or infall region of the cluster, that cannot be 
explained by our ram-pressure stripping model, or the ef¬ 
fect of backsplashing. This finding suggests that an addi¬ 
tional mechanism must also be at play. One possibility is 
that galaxies are pre-processed in their previous enviroment 
(e.g. smaller groups, or the field) before they fall into clus¬ 
ters. We will investigate this in more detail in future papers. 
We conclude, however, that if a blue, gas-rich, star-forming 
galaxy has not been quenched by the time it falls into the 
cluster, it is most likely to experience gas removal as it first 
passes through the cluster’s ICM due to the effect of ram- 
pressure. 

Finally, this work highlights the usefulness of phase- 
space diagrams to infer the orbital histories of galaxies in 
regular clusters, and the effect of cluster processes such as 
ram-pressure stripping. It should be noted, however, that 
this method cannot be robustly applied to very irregular 
(e.g. merging) clusters. Moreover, our phase-space analy¬ 
sis proved to be particularly insightful when using sensitive 
tracers of environmental effects such as HI, in combination 
with optical and UV data. 
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APPENDIX A: SENSITIVITY OF THE 
RAM-PRESSURE STRIPPING MODEL IN PPS 

In this section we study the sensitivity of our ram-pressure 
stripping model (Section 14.Ill with variations of the cluster 
and galaxy model parameters. 

We take A963_l and start by considering the case of a 
MW galaxy that is completely stripped: truncation radius 
r' = 0 and 77 = 1 , which corresponds to the dashed green 
line in the top panel of Figure [Xl] (same as dashed grey line 
in the left panel of Figure Q • This line separates the region 
where P„ ra is able to strip all the gas in the modeled galaxy 
(left side of the line, P ra m > n sa ;) from the area where the 
anchoring force of the galaxy is still able to hold on to some 
of the gas (right side, P ra m < n ga i). 

Now, what happens if P ra m is bigger or smaller than 
the value computed for A963_l? Or when H ga i differs from 
that of the MW? We illustrate this with the different lines 
in the top panel of Figure [ATI The lines correspond to dif¬ 
ferent values of 77, as indicated by the legend. For example, 
the blue line shows the region where Pram is 10 times higher 
than in A963_l (alternatively, it could also be the case of a 
galaxy with H ga i 10 times smaller than the MW, in a A963_l 
cluster). The red line instead corresponds to a cluster with 
ram-pressure 5 times weaker. This set of values attempts to 
account for the variation of U ga i and p ICM . As the figure 
shows there is a very small amount of variation, when con¬ 
sidering variations of the order of 10 (blue and red lines). 
Significant deviations from the model only arise in very ex¬ 
treme cases (i.e. when 77 is modified by a factor of 100 ). 
We thus conclude that the region where ram-pressure is ex¬ 
pected to dominate (dashed greeen area in Figure [7]) will 
not change significantly if the parameters used to model our 
cluster or our test galaxy are somewhat different. 

Because not all of our galaxies are like the MW, we 
specifically test the sensitivity of our model (again, for 
A963_l) with galaxies of different surface density, which de¬ 
pends on R, 1, and M*, as shown by equation 1101 

We take ave r age v alues of Rd in 3 bins of stellar mass 
from lFathi et al] 1 20 id . see Figure 15) and compute n ga i in 
each case. The resulting ram-pressure stripping model (for 
77 = 1 ) is shown by the thick solid lines on the bottom panel 
of Figure lATl Notice the stellar mass range considered is 
very large, from < 1 O 8 M 0 , to 10 11 — 10 12 Mq. In this plot, 
our MW model (M* = 4.6 x 1O 1 O M 0 ) would lie somewhere 
between the blue and red line. We further accounted for the 
scatter or uncertainty in Rd in each case, by computing the 
associated uncertainty in each model (In, shaded regions). 
The values used in each model are: 

• For galaxies with M* < 1O 8 M 0 , the average scalelength 
is (Rd) = 0.238±0.094 kpc (green) 

• For M* = 10 9 — 1 O 1 O M 0 , (Rd) = 1.52±0.65 kpc (blue) 

• For M* = 10 11 - 1O 12 M 0 , (R d ) = 5.73± 1.94 kpc (red). 
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Ram-pressure model: rj variations 



r /^200 


Stripping in different galaxies 



Figure Al. Top: The ram-pressure stripping model for a MW 
galaxy in A963-1 is shown for different ram-pressure stripping in¬ 
tensities (77 = Pram/^lgal)} as labeled. Bottom: The ram-pressure 
stripping model for cluster A963_l (using 77 = 1 ) is tested against 
galaxies of different stellar mass, as labeled. The thick lines corre¬ 
spond to the model using the average scale length for that mass 
range (see text), and the shaded areas show the uncertainties 
introduced by the scatter in Rd- In both panels the axis range 
mimics that of Figure [7] (left panel) for easy comparison. 


T hese values were t aken directly from the text in Section 4.3 
of lFathi et alj (l20ld ). We did not use different mass bins to 
avoid interpolating (Rd) from their figure. The plot shows 
that for galaxies more mass ive than 10 10 M(7), and Rd values 
within lcr of those quoted in lFathi et alj J20ld ). the stripped 
region PPS is very similar. 

We thus conclude that, there is little variability induced 
by different cluster models, and that for a mass-limited sam¬ 
ple, the method is very robust in predicting the location in 
PPS where ram-pressure is expected to strip off the gas in 
galaxies. 
















